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2 1- Introduction 



ABSTRACT 



X-ray variability of active galactic nuclei (AGN) and black hole binaries can be analysed by means of the power spectral density 
(PSD). The break observed in the power spectrum defines a characteristic variability timescale of the accreting system. The empirical 
variability scaling that relates characteristic timescale, black hole mass, and accretion rate (Tg oc Ml^/A/"'*) extends from supermas- 
sive black holes in AGN down to stellar-mass black holes in binary systems. We suggest that the PSD break timescale is associated 
with the cooling timescale of electrons in the Comptonisation process at the origin of the observed hard X-ray emission. We find that 
the Compton cooling timescale directly leads to the observational scaling, which naturally reproduces the functional dependence on 
black hole mass and accretion rate (tc »; Mj^^/M). This result simply arises from general properties of the emission mechanism and 
is independent of the details of any specific accretion model. 
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, X-ray variability is a characteristic property of all accreting 
black hole systems, including supermassive black holes in active 
galactic nuclei (AGN) and stellar-mass black holes in binary sys- 
, tems (black hole binaries, BHB). The observation of rapid and 
' aperiodic variations indicates that X-ray emission is produced 
in the innermost regions of the accretion flow, close to the cen- 
tral black hole. The study of X-ray variability is thus thought to 
provide information on the physical emission mechanisms and 
I possibly on flie nature of the underlying accretion flows. 

Quantitative analysis of X-ray variability is often performed 
by means of the power spectral density (PSD). The observed 
. variability power spectrum can be roughly modelled by power 
■ laws of the form Py oc v", where v - l/T is the temporal fre- 
quency and a the power law slope. A typical AGN power spec- 
trum is characterised by a slope of a ~ -2 at high frequen- 
cies (corresponding to short timescales), flattening to a slope of 
a ~ -1 at lower frequencies (or equivalently longer timescales). 
A similar shape is observed in the power spectra of galactic 
black holes in the high/soft state. In addition, a second break 
to a slope of a ~ 0, occurring at still lower frequencies, is seen 
in the low/hard state of black hole binaries. In the following we 
consider the high-frequency PSD break, vb (from a ~ -1 to 
a 2), in both AGN and black hole binaries. The correspond- 
ing break timescale, Tg = I /vg, is considered to be a character- 
istic X-ray variability timescale of the accreting system. Typical 
break timescales measured in AGN are of ~day-timescales, 
while the corresponding break timescales seen in black hole bi- 
naries are of the order of fractions of seconds (1VIcHardvl l2010t 
[Gilfanov 2010). The difference in break timescales observed be- 
tween supermassive and stellar-mass black holes is broadly con- 
sistent with a simple scaling with black hole mass, and linear 



mass-timescale relatio ns have been proposed (iMarkowitz et al.l 
l2003l:IPapadakisll200l . 

A more quanti tative X -ray variability scaling has been ob- 
tained by McHard v et al.l (|2006), who considered the explicit 
dependence of the characteiistic timescale on both black hole 
mass and accretion rate. A major result is that this observa- 
tional scaling relation extends from supermassive black holes in 
AGN down to stellar-mass black holes in binary systems, cover- 
ing many orders of magnitude in black hole mass. The fact that 
both AGN and black hole binaries seem to follow the same em- 
pirical scaling has been claimed to be a strong support for the 
scale invariance of accreting black hole systems, and led to the 
suggestion that AGN are simply scaled-up galactic black holes 
(McHardv et al. 2006). This in turn led to the proposal of the 
so-called 'vaiiability fundamental plane', which relates charac- 
teristic timescale, black hole mass, and accretion rate. 

However, the physical oiigin of the PSD break is not defini- 
tively established yet. Vaiious inteipretations of the PSD break 
timescale have been proposed in the literature, trying to link 
the observed break timescales with different accretion model 
timescales. At first, the measured break timescales have been 
compared with characteristic physical timescales of the accre- 
tion disc, inclu ding orbital, thermal, and viscous times cales at 
different radii ('Marko witz et al.ll2003l lPapadakisir2004 . More 
recently, the PSD break has been discussed in the framework 
of the inner propagating accretion flow fluctuations model, in 
which fluctuations originating at different outer radii propa- 
gate inwards, leading to modulations of the accretion rate at 
inner radii and hence to variations in the observed X-ray flux 
(Lvubarskii 1997; Churazov et al. 2001). The optically thick ac- 
cretion disc is assumed to be truncated at some distance from 
the centre and the inner region is replaced by a hot, optically thin 
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flow component. In this picture, the break timescale is associated 
with a characteristic timescale, usually the viscous time, at the 
truncation radius of the standard accretion disc. In an alternative 
approach, we have discussed the characteristic X-ray variability 
timesc ale within the framework of a cl umpy accretion model in 
AGN dlshibashi & Courvoisieij|2009bl) . In this scenario, accre- 
tion onto the central black hole occurs through a sequence or 
cascade of shocks, which are at the origin of the observed radi- 
ation. Accreting elements first interact in the form of optically 
thick shocks, giving rise to optical/UV emission, later followed 
by optically thin shocks giving rise to X-rays through inverse 
Compton emission. Optical/UV photons emitted in the optically 
thick shocks form the seed photons for Compton-cooling of the 
electrons heated in the optically thin shocks. The characteristic 
heating/cooling timescale of the electrons in the optically thin 
shocks has been associated with the PSD break timescale. 

Although the exact geometry of the inner coronal flow 
is still unknown, hard X-ray emission is widely attributed to 
Compton processes. In the inverse Compton process, soft seed 
photons (possibly emitted as thermal emission from the accre- 
tion disc) are upscattered to higher energies by energetic elec- 
trons in a hot, optically thin corona. The observed hard X- 
ray power law component of accreting black holes is indeed 
well represented by Comptonisation models (iTitarchukI 1 1 994L 
and references therein). As hard X-rays are produced by the 
Comptonisation mechanism, it seems natural to consider physi- 
cal timescales associated with the emission process itself as pos- 
sible candidates for the identification with the observed X-ray 
PSD break timescale. In the following we consider the charac- 
teristic timescale associated with the Compton cooling of the hot 
electrons. 



2. Compton cooling timescale 

The power emitted by a single electron in the inv erse Compton 
scattering is given by (iRvbicki & Lightmanll 19791) 



^c = TCo-T7^/3^Uy , 



(1) 



where (Tt is the Thomson cross section and Uy the photon energy 
density. 

The Compton cooling timescale, which sets the character- 
istic timescale for the cooling of electrons by inverse Compton 
emission, is defined as 



where ri~0. 1 is the canonical radiative efficiency. The radial dis- 
tance R can be expressed in units of the Schwarzschild radius as 
R - ^Rs , where Rs - 



ICMb 



The Compton cooling timescale (Eq. O is now written as 



fc = 



a-jc'^ J] M 



(5) 



We see that the Compton cooling time is mainly determined 
by two of the black hole fundamental parameters, i.e. mass and 
accretion rate. As the Compton cooling time is directly linked to 
the physical process of X-ray emission, it can be considered to 
be a characteristic X-ray timescale of the accreting system. 

For a typical AGN of mass Mbh ~ IO^Mq and accreting in 
the radiatively efficient regime (with typical Eddington ratios of 
L/Le ~ 0.1), the numerical value of the characteristic timescale 
(Eq.|5]l is of the order of 



fAc^~0.l77„',4„M:]M2d, 



(6) 



where Mbh = M? • IO^Mq and M = M_2 ■ O.OlMo/yr. We 
see that the characteristic Compton cooling timescale in AGN 
is of ~day-timescales. We can also estimate the Compton cool- 
ing timescale in black hole binaries by introducing parameters 
appropriate for galactic black holes in Eq. (|5). Considering an 
object of about Mbh ~ IOMq accreting at a few percent of the 
Eddington limit, and taking into account the fact that the seed 
photon luminosity is only a fraction of the bolometric luminos- 
ity, the characteristic timescale is of the order of Ibhb ~ 0.07s. 

We observe that the Compton cooling timescale is shorter for 
smaU black hole mass and/or high accretion rate, scaling as 



tc oc 



Kh 
M 



(7) 



The Compton cooling of the hot electrons leads to X-ray 
emission. The average electron energy can be estimated by 
assuming an equilibrium situation between Compton coolin: 



i 



and Coulomb heating rates: Lc = Lh (Is hibashi & Courvoisie: 
I2009ah . The resulting electron energy scales as Eg oc nj Uy , 
where iie is the electron number density. The average X-ray lu- 
minosity is given by the Compton luminosity of a single electron 
(Eq.[l]l multiplied by the average number of electrons: 



tc 



Ee 



(2) 



{Lx) ~Lc-Ne^ nl^'uf , 



(8) 



where E^ is the average electron energy. In the non-relativistic 
limit (y ~ 1, y6^ oc Ee), the dependence on the electron energy 
cancels, and Eq. (|2| leads to 



tc 



3meC 1 

SCTj Uy 

The radiation density at distance R is given by 



AjtR^c 



(3) 



(4) 



where L, is the seed photon luminosity for inverse Compton 
scattering. Assuming that the seed photons are generated in the 
radiatively efficient accretion process, the seed photon luminos- 
ity is directly proportional to the accretion rate, L, - rjMc^, 



and we see that the dependence on the electron number den- 
sity is stronger than that on the photon energy density. Thus 
X-ray luminosity variations are mainly driven by electron den- 
sity modulations associated with the inhomogeneities of the up- 
scattering medium. In the clumpy accretion scenario, the origin 
of the density modulations was associated with the inhomoge- 
neous structure of the accretion flow resulting fr om the succes- 
sion of shocks (ishibas hi & Cour voisier"2009b'). Here we note 
that the physical origin of the modulations can be diverse and 
does not need to be specified. Density fluctuations on timescales 
shorter than the electron cooling time (tf < tc), which at equi- 
librium is equal to the electron heating time, cannot be propa- 
gated and are smoothed out. The suppression of fluctuations on 
timescales shorter than a critical value, given by the Compton 
cooling time, induces a break in the observed variability power 
spectrum, which may be associated with the PSD break. 
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3. Comparison with tlie PSD brealc timescale 

Observational measurements of the high-frequency PSD break 
(from a ~ -1 to o- ~ -2) have been performed for 
a number of AGN and galactic black holes. The observed 
break timescales in Se yfert galaxies are usua l ly in the range 
0.01 to a few davs (JMarkowitz eTaP 120031: iMcHardv etaP 



l2004t lUttlev & McHardy||2005h . while the corresponding break 
timescales measured in black hole binaries are of the order of 
seconds or less (Gilfanov 2010, and references therein). The 
Compton cooling timescales estimated in the previous section 
(~day-timescales in AGN and fractions of seconds in BHB) 
thus correspond to the typically measured values of the PSD 
break in both AGN and galactic black holes. The character- 
istic timescale therefore seems to scale roughly linearly with 
black hole mass, from stellar-mass black holes up to supermas- 
sive black holes. However, an additional dependence on the ac- 
cretion rate is required by more detailed observations of AGN: 
for a given black hole mass, higher accretion rate objects have 
shorter break timescales t han lower accretion rate counterparts 
(lUttlev & McHardvll2005h . The empirical relationship between 
variability timescale, black hole mass, and accretion rate has 
been obtained by McHardv et al. (2006) based on a main sample 
of AGN including two galactic black holes in radio-quiet states. 
The resulting empirical scaling can be written as 



M 



1.12 



1.0.98 



(9) 



where niE - Lboi/LE is the ratio of the bolometric luminosity 
to Eddington luminosity. This observational scaling is shown to 
extend from AGN down to stellar-mass black holes in binary 
systems. As the Eddington luminosity is given by the black hole 
mass (Le oc Mbh) and assuming that the bolometric luminosity 
is proportional to the accretion rate (Lhoi oc M), the empirical 
scaling (Eq.|9]) can be re-expressed directly in terms of the phys- 
ical parameters as 



r„~ 



MIL 

M0.98 



(10) 



M-„u 



Comparison of the Compton cooling timescale, tc o^ —ff- 
(derived in Eq.|5]l, with the observational scaling (Eq.fTOli shows 
a very close match. We observe that the characteristic timescale 
associated with the Compton cooling process directly repro- 
duces the observed dependence on black hole mass and accretion 
rate. We note that a similar result has been previously derived but 
within the specific framework of the c lumpy accretion model in 
AGN (llshibashi & Courvoisiedl2009bl) . 



4. Spectral-timing relation 

X-ray variability properties are also closely coupled with X- 
ray spectral properties. Observations of 'spectral-timing' corre- 
lations have been reported in both AGN and black hole binaries, 
whereby an increase in the characteristic frequency is correlated 
with a steepening of the X-ray spectrum. The characteristic noise 
frequency is observed to increa se with the steepening of the en- 
ergy spectrum in Cygnus X-1 (iGilfanov et al.lll999l) . and sim- 
ilar patterns are also observed in other X-ray binary systems, 
e.g. GX 339-4 (Revnivtsevet al] l200lh . In the case of AGN, 
the average X-ray spectral slope is found to be positively cor- 
related with the normalized characteristic frequency defined as 



Vnorm = Vfi X M bh, where vb is the br eak frequency observed in 
the X-ray PSD (Paoada kis et"ani2009h . 

It is known that Comptonisation processes lead to power law 
spectra, and the resulting X-ray spectra observed in the hard en- 
ergy band are usually described by power laws, characterised by 
the photon index F. A significant correlation between photon in- 
dex and Eddington ratio has been reported in numerous AGN 
samples (Shemmeret al. 2008; Kelly et al. 2008; Risal iti et all 
2009), suggesting that the X-ray spectral slope is correlated with 
some form of accretion rate. Within the Comptonisation sce- 
nario, the physical interpretation of the observed T-L/Le corre- 
lation is that an increase in the accretion rate leads to an enhance- 
ment in the seed photon emission, which in turn leads to effi- 
cient Compton cooling of the hot corona, which finally leads to 
a steepening of the hard X-ray spectrum. On the other hand, we 
have seen that the Compton cooling time scales inversely with 
the seed photon energy density (Eq. [3]l, thus an increase in the 
seed photon emission also leads to a decrease in the characteris- 
tic timescale. Combining the two results we may account for the 
observed correlation between characteristic frequency and spec- 
tral steepening. Therefore in the context of Comptonisation, both 
the X-ray characteristic timescale and the spectral slope seem to 
be determined by the supply of seed photons, which in turn de- 
pends on the accretion rate. Hence, for a given black hole mass, 
higher accretion rate systems tend to have both steep spectra and 
short characteristic timescale s, in agreement with th e observed 
'spectral-timing correlation' (iPapadakis et al.ll2009i) . This sug- 
gests that both timing and spectral properties of accreting black 
hole systems are mainly determined by the accretion rate param- 
eter. 

5. Discussion 

The break observed in the X-ray PSD is thought to provide 
insight into the underlying variability mechanism and various 
physical interpretations of the break timescale have been pro- 
posed. As a first approach, the observed PSD break timescales 
have been compared with physical timescales of the accretion 
disc in the innermost regions of the accretion flow where X- 
ray emission is thought to be produced. The shortest timescale 
is given by the dynamical or orbital timescale, t^. This is of 
the order of ~hours for typical AGN parameters and thus too 
short compared to the break timescales of the order of ~days 
observed in Seyfert galaxies. On the other hand, for a stan- 
dard geometrically thin disc (with aspect ratio H/R « 1) the 
viscous timescales (of the order of >years) are found to be 
a few orders of magnitude lo nger compared to the observed 
values (Markowit z et al.ir2003h . The closest values to the mea- 
sured break timescales are perhaps given by thermal t imescales, 
and t hermal instabilities have also been considered (iPapadakisI 
12004 *). We note that the physical timescales of the accretion disc 
are all proportional to the black hole mass and thus the approxi- 
mate scaling of the characteristic timescale with black hole mass 
is easily obtained; however, the dependence on the accretion rate 
is less obvious and cannot be simply accounted for 

The currently favoured interpretation of the PSD break is 
based on the inner propaga ting fluctuations model (Lvubarsldj 
1 19971 : IChurazov et al.ll2()()ll) . The variability power is found to 
decrease on short timescales, with a rapid fall-off below the 
break, suggesting that the break timescale might be associated 
with some cut-off in the accreting system. Accordingly, the 
break timescale has been associated with the viscous timescale 
(tv ~ mass (H/R)^) at the inner edge of the optically thick ac- 
cretion disc. The accretion disc is assumed to be truncated at 
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some distance from the centre, the location of the truncation ra- 
dius being connected with the accretion rate: at high accretion 
rates the optically thick disc extends down to the last stable or- 
bit, while at low accretion rates the d isc is truncated at large 
distances (IZdziarski & Gierliiiskill2004l and references therein). 
In this picture, the value of the characteristic timescale is related 
to the position of the disc truncation radius. If the disc is trun- 
cated far away from the centre, the associated viscous timescale 
is long; while with increasing accretion rate, the truncation ra- 
dius becomes smaller and the associated characteristic timescale 
correspondingly shorter Once the accretion disc has reached the 
last stable orbit, a further decrease of the viscous timescale can 
only be obtained by varying the disc thickness parameter (H/R), 
assu ming that the H/ R ratio increases with increasing accretion 
rate (lMcHardvll20 1 0*) . Therefore by combining variations of the 
disc truncation radius with variations of the disc thickness pa- 
rameter one may qualitatively account for the observed trend 
with accretion rate. 

Although the PSD break timescale has usually been associ- 
ated with some characteristic timescale of the accretion disc, it 
is well known that the accretion disc itself cannot account for 
the hard X-ray emission component of accreting black holes. 
Consequently, the identification of a characteristic disc timescale 
at a particular radius with the characteristic X-ray variability 
timescale is not straightforward, and other interpretations may 
be considered. On the other hand, it is widely agreed that hard X- 
ray emission is produced by Comptonisation processes. It seems 
therefore interesting to see whether there is any characteristic 
timescale associated with the emission process itself that might 
be compared with the observed PSD break timescale. 

We have previously associated the PSD break timescale 
with the heating/cooling time of the electrons in the op- 
tically thin shocks w ithin the clumpy accretion scenario 
(llshibashi & Courvoisi er 2009b). Here we suggest that the X- 
ray PSD break timescale could be associated with the cooling 
time of electrons in the Comptonisation process, independently 
of the details of any specific accretion model. The characteris- 
tic X-ray variability timescale is simply given by the Compton 
cooling time; the observed scaling with black hole mass and 
accretion rate is quantitatively reproduced, without additional 
assumptions. It is quite remarkable that the observational scal- 
ing simply derives from the definition of the Compton cooling 
timescale with a basic assumption of a radiatively efficient ac- 
cretion and a quasi-spherical distribution of the seed photon pop- 
ulation. 

The Comptonisation process is thought to operate in a 
similar way in s upermassive and stellar-mass black holes 
(iGliozzi et alJl2011 ). We have seen that numerical estimates of 
the Compton cooling time give correct order-of-magnitude val- 
ues of the PSD break timescales in AGN and binary systems. In 
our picture, the characteristic timescale is set by the Compton 
cooling time in both AGN and galactic black holes, suggest- 
ing that the underlying emission mechanism is scale-invariant. 
In addition, if the characteristic X-ray variability timescale is 
given by the Compton cooling timescale, this may additionally 
support the hypothesis that the main mechanism of hard X-ray 
emission is indeed Comptonisation in both classes of accreting 
black holes. 



tion disc, such as the viscous timescale at the truncation ra- 
dius. Here we suggest another possibility and associate the break 
timescale with a physical timescale directly linked to the X-ray 
emission mechanism, i.e. the cooling time of the electrons in 
the Comptonisation process. The association seems natural since 
hard X-rays are generated in the Comptonisation process and it 
also provides a possible physical interpretation of the PSD break. 
We obtain that the Compton cooling time directly reproduces the 
correct functional dependence on black hole ma ss and accretion 
rate, d efining the empirical variability scaling of iMcHardv et al.l 
(l2006h . We show that this result is model-independent while the 
detailed physical realisations of the accretion flows might take 
different forms in different accreting systems. 
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6. Conclusion 

Various physical interpretations of the PSD break have been pro- 
posed in the literature. The observed break timescale is usu- 
ally associated with some characteristic timescale of the accre- 



